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Abstract. Using the Riemann-Hilbert approach, the ^-function correspond- 
ing to the solution of the first Painleve equation y xx = 6y 2 + x with the 
asymptotic behavior y ~ ±v _ x/& as \x\ — > oo is constructed. The exponen- 
tially small jump in the dominant solution and the coefficient asymptotics in 
the power-like expansion to the latter are found. 



1. Introduction 
The Painleve first equation 1 

(-Pi) Vxx = 6y 2 + x, 

is the simplest of the six classical equations of Painleve-Gambier and can be 
derived from any other Painleve equation using certain scaling reductions 3 . The 
recent interest to this equation is due to its significant role in various physical 
models. 

For instance, equation Pi describes certain solutions to KdV and Bussincsq equa- 
tions @I E], bifurcations in some non-integrable nonlinear models 6 , continuous 
limits in matrix models of quantum gravity 7, ,8, S| HOj ; ^-function associated with 
Pi appears in n-large asymptotics of semi-classical orthogonal and bi-orthogonal 
polynomials jlllll2j and thus becomes a primary object in the problem of Laplacian 
growth |13| . 

In context of the string theory, "physical" solutions of Pi are distinguished 
from "non-physical" ones by the monotonic asymptotic behavior as i — * — oo 
|U 1141 115| . There are two kinds of such monotonic boundary conditions, 
i.e. y(x) ~ x/6. Using elementary perturbation analysis, the solution y(x) = 

is unique as being a background to a 2-parametric family of 
oscillating solutions. Solutions approaching a positive branch of the square root as 
x — > — oo, i.e. y{x) ~ y/—x/6, form a 1-parametric family parameterized by an am- 
plitude of an exponentially small perturbation to a power-like dominant solution. 
These solutions have the asymptotic expansion y(x) = ^—x/6 J2kLo a k{—x)~ 5k / 2 + 
0(x~°°), whose coefficients ak, admit a combinatorial interpretation [UAH?] . 

In the problem of Laplacian growth without surface tension (Hcle-Shaw problem, 
quantum Hall effect etc.), the shape of a growing droplet is described using certain 
^/-function, see [T^J. If the droplet develops a cusp singularity, this ^/-function can 
be approximated by a ^/-function associated to the first Painleve transcendent, [T%| . 
Certain "physical" asymptotic conditions imposed on this ^-function determine 
the relevant Stokes multipliers Sk- In turn, these Sk pinch out the monotonic as 
x — > — oo Painleve function y(x) ~ y/—x/6. 
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Equation Pi has unexpectedly rich asymptotic properties in the complex x- 
plane. P. Boutroux ^H] has shown that, generically, asymptotics of the Painleve 
first transcendent as \x\ — ► oo is described by the modulated WeierstraB elliptic 
function whose module is a transcendent function of arg x. Furthermore, the module 
function is such that the elliptic asymptotic ansatz degenerates along the directions 
arg a; = ir+^-n, n = 0, ±1, ±2. P. Boutroux called the corresponding trigonometric 
expansions "truncated" solutions. Their 1- and 0-parameter reductions, if they 
admit analytic continuation into one or two neighboring sectors of the complex 
s-plane, were called by Boutroux "bi-truncated" and "tri-truncated" solutions. All 
bi- and tri-truncated solutions have the algebraic leading order behavior, y(x) ~ 
±y/—x/6, perturbed by exponential terms. 

We call a discontinuity in the asymptotic form of an analytic function the Stokes 
phenomenon. In the case of Pi, a jump in the phase shift of a modulated elliptic 
ansatz across the rays arg x = tt + ^-n is called the nonlinear Stokes phenomenon. 
For bi- and tri-truncated solutions, a jump in the exponentially small perturbation 
of a dominant solution resembles the well known Stokes phenomenon in the linear 
theory and thus is called the quasi-linear Stokes phenomenon. 

In |201 1211 122j . equation Pi was studied further using classical tools like the 
perturbation approach and the method of nonlinear integral equations. Mainly, 
these articles discuss the behavior of the Painleve transcendents on the real line. 
The recent paper |23| adopts the same approach carefully studying the behavior of 
the tri-truncated solution on the negative part of the real line. 

In |241 125| . the multiple scale analysis was applied to Pi (and P?) to find a 
precise form of the phase shift in the elliptic asymptotic ansatz within complex 
sectors between the indicated rays. In |14U15j . the Witham averaging method was 
used to describe the elliptic tail of the monotonic at — oo solution of Pi. 

The isomonodromy deformation approach to Painleve equations, see |2(j| 1271 125] . 
was applied to equation Pi in |291 131)1 131) . In this way, the asymptotics of the 
Painleve functions is expressed in terms of the Stokes multipliers of an associated 
linear system. Then the equation of a monodromy surface yields connection formu- 
las for the asymptotic parameters along different directions of the complex x-plane. 
A complete description of the nonlinear Stokes phenomenon in Pi is given in |3()j . 
A heuristic description of the quasi-linear Stokes phenomenon in Pi can be found 
in E!]. 

Using the Borel transform technique and some assumptions on the analytic prop- 
erties of the relevant Borel transforms, as well as the isomonodromy deformation 
approach based on the so-called exact WKB analysis, Y. Takei (22] has re-derived 
the latter result (look for more discussion in |33|1. 

In the present paper, we construct the ^-function associated with the monotonic 
as \x\ — ► oo solution of Pi and rigorously describe the relevant quasi- linear Stokes 
phenomenon. Our main tool is the Riemann-Hilbert (RH) problem. We observe 
that the jump graph for our RH problem can be decomposed into a disjoint union of 
two branches, one of which is responsible for the background y^—x/6 while another 
one produces the exponentially small perturbation of the dominant solution (look 
for similar observation in P2 case). Using the steepest descent approach of 
Deift and Zhou |34|. we prove the unique solubility of this problem and compute 
the asymptotics of the Painleve transcendent. 
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Applying a rotational symmetry, we prove the existence of five solutions, yi n (x), 
n = 0,±1,±2, asymptotic to \J e~ i7r x/6 as \x\ — > oo in the respective overlapping 
sectors argrr € (-f - 4fn, ?f - ^fn), see ^IW^ . |2"7r2")) . and find the exponentially 
small differences y4( n -i){x) — U4n(x), n — 0, ±f , ±2, see (|2.7fll . The latter constitute 
the quasi-linear Stokes phenomenon. 

A collection of the functions y4 n (x), n — 0,±1,±2, forms a piece-wise mero- 
morphic function y(x) ~ \J e~ in x/6 as — * oo. The moments of this function 
immediately yield the asymptotics as k — > oo for the coefficients a k l|3.12|l of the 
x-series expansion to the dominant solution (|3 . 1 \ . 

For the first time, the formula for the coefficient asymptotics was found in |32| . 
The authors of [23] studied the recurrence relations for the same coefficients by 
direct means and prove a similar asymptotic formula modulo a common factor (an 
advanced version of the direct approach to a generalized recurrence relation which 
contains one for Pi as a special case can be found in [HE])- The exact value of this 
common factor was announced in [2J with the reference to the method of (36[ based 
on the Borel transform formula. In contrast, we do not use the Borel transform 
technique at any stage of our investigation. 

The paper is organized as follows. In Section |2 we recall the Lax pair for Pi, 
formulate the relevant RH problem and solve it asymptotically. Using the approx- 
imate ^-function, we find the asymptotics of the bi- and tri-truncated Painleve 
transcendents and of the relevant Hamiltonian functions. In Section UJ we find the 
coefficient asymptotics in the power- like expansion to the formal solution of Pi . 



2. RlEMANN-HlLBERT PROBLEM FOR PI 

Introduce generators of sm(2, C), 03 — a + = (q 1 ), er_ = (i°) an d 

the Pauli matrices oi = a + + cr_ and 02 = + iff- together and consider the 

system of matrix differential equations for 'J, see [371 138j . 

(2.1a) =A(X,x) = -za 3 + (2A 2 + 2y\ + x + 2y 2 )a+ + 2(A - y)cr_, 

CM 

(2.1b) ^U" 1 =U{\,x) = {\ + 2y)o + + o_. 
ox 



Compatibility of (|2. lafl and (|2.1b|l implies that the coefficients z and y depend on 
the deformation parameter x in accord with the nonlinear differential system 



(2.2) 



Vx = z, 

z x = 6y 2 + x, 



which is equivalent to the classical Painleve first equation |Pi| Following 1201, see 
also |39j , linear equation (|2.1al) has the only one (irregular) singular point at infinity, 
and there exist solutions 'J'fc(A) of H2.1|l with the asymptotics 
(2.3) 

tf fc (A) = \^^{^ + ai){l-'K<j ? ,\- 1 ' 2 + 0{\- 1 ))e e ^>\ 9{\) = |A 5 / 2 + XX 1 / 2 , 
as 

(2.4) A^oo, Xen k = {XeC: argAe (f (fc - |), ^(fc + |))}, k E Z. 

Solutions ^fe(A), k G Z, are called the canonical solutions, while sectors flk are 
called the canonical sectors. Canonical solutions ^(A) are uniquely determined by 
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(|2.3(l - (|2.4(l and solve both equations l|2.1J) . They differ from each other in constant 
right matrix multipliers Sk called the Stokes matrices, 

(2.5) *fc +1 (A)=* fe (A)S fc , Si*-i=(j S2fc 1 " 1 ), S ™=( s ] k o 

Observing that all solutions of p.lafl are entire functions, thus 

(2.6) * fc (e 2 "A) = * fc (A), 
and using the relation 

(2.7) *A:+5(e 27rl A) = ty k (\)iai 

which follows from the definition of the canonical solutions and the asymptotics 
(|2.3I) . 1)2. 4|l . we readily find the constraints for the Stokes matrices |29) . 

(2.8) = CTiSfcCTi, S , iS , 2 S , 35'45'5 = ici, 
or, in the scalar form, 

I) Sfe+5=Sfc, 1 + SfcSfc+i = -isk+3, keZ. 



Thus, generically, two of the Stokes multipliers Sk, k G Z, determine all others. 

The inverse monodromy problem consists of reconstruction of ^fe(A) using the 
known values of the Stokes multipliers Sk- It can be equivalently formulated as 
a Riemann-Hilbert (RH) problem. With this aim, introduce the union of rays 
7 = P U (U| =1 7fc_ 3 ) , where lk = {A e C : arg A = ^*}, fc = -2, -1, 0, 1, 2, and 
p = {A € C: arg A = 7r}, all oriented toward infinity. Denote the sectors between 
the rays p and 7_ 2 by uj_ 2 , between jk-i and 7/., fc = —1,0,1,2, by w^, and 
between 72 and p by W3. All the sectors are in one-to-one correspondence with 
the canonical sectors 1)2. 4f> . see Figure ^ 

Let each of the sectors u>k, k = —2,-1,. ..,3, be a domain for a holomorphic 
2x2 matrix function ^/^(A). Denote the collection of 4 r fc(A) by ^(A), 

( 2 - 9 ) *( A )Le„ fc = ^( A )- 

Let \E f _|_ (A) and \&-(A) be the limits of ^E'(A) on 7 to the left and to the right, 
respectively. 

Let 0(A) = I A 5 / 2 + xX 1 / 2 be defined on the complex A-plane cut along the 
negative part of the real axis. The RH problem we talk about is the following one: 

(1) Find a piece- wise holomorphic 2x2 matrix function \& (A) such that 

(2.10) lim A 1/2 (-L(cr 3 + <7 1 )A~* <T3 *(A)e- eo ' 3 - i) exists and is diagonal; 

(2) on the contour 7, the jump condition holds 

(2.11) tf + (A) = ¥_(A)S(A), 

where the piece-wise constant matrix S{\) is given by equations: 

(2.12a) S(X)\ = S k , S 2 k-1 =I + s 2 k-l0-+, S 2 k = I + s 2 kcr~, 

(2.12b) S(A)| p = -**!, 

with the constants s k satisfying the constraints (JS^Sfll; 
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FIGURE 1 . The RH problem graph for Pi 



oo, 



Because VP(A) satisfies the asymptotic condition 

(2.13) Y(X) := -±(a 3 + ^JA" W(A)e- flos = 

= (i-^ + ^ + o(x-^) ^ + o(x~^) \ 
\ ^ + o(x-^) 1 + ^ + ^l + o(x-^))' 

where 

(2.14) J£ = \z 2 - 2y 3 ~ xy, 

the solution y(x) of the Painleve equation can be found from the "residue" of Y(X) 
at infinity, 

(2.15) y = 2 lim AYi 2 (A) = 2 lim AF 2 i(A). 

A — >oo A — >oo 

Remark 2.1. It is easy to see that IH is nothing but the Hamiltonian for the Painleve 
first equation with the canonical variables q = y and p = z. 

Equation 12.15fl specifies the Painleve transcendent as a function y = f(x, {sk}) 
of the deformation parameter x and of the Stokes multipliers Sk ■ Using the solution 
y = f(x,{sk}) and the symmetries of the Stokes multipliers described in |29| . we 
obtain further solutions of|Fi| 



(2.16a) y = f(x, {-«=]?}), 

(2.16b) y = e^ n f(e l ^ n x,{s k+2n }) 7 n g Z, 

where the bar means the complex conjugation. 
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By technical reason, to find the asymptotics of y(x), we use below not Y(X) 
but related auxiliary functions %(A) and X(X) with expansions (|2.45|) and Ij2.50|l . 
respectively. The latter involve differences y — y(x), y(x) are known, which can be 
estimated using singular integral equations with contracting operators. 

2.1. Asymptotic solution for sq = 0. Let us consider the RH problem above 
where so = assuming that \x\ — > oo within the sector arg 2 £ [%"> 7r ]- Equations 
(|2.12a|l imply that "J (A) has no jump across the ray 70 = {A £ C: arg A = 0}. The 
constraints ij^I)) reduce to the following system of equations, 

(2.17) s_ 2 = s 2 = s_i +si=i. 

Our first step in the RH problem analysis consists of introduction an auxiliary 
(?-function, 

(2.18) 5 (A) = !(A + 2A ) 5 / 2 -4A (A + 2A ) 3 / 2 , A = y/e~^x/6, 

defined on the complex A-plane cut along the ray (—00, — 2Ao]. The asymptotics of 
g-function as A — > 00 coincides with the canonical one, 

(2.19) g(\) = |A 5 / 2 -6A§A 1 / 2 -4AgA- 1 / 2 + 0(A- 3 / 2 ) = |A 5 / 2 + ^A 1 / 2 +0(A- 1 / 2 ). 

Let us formulate an equivalent RH problem for the piece-wise holomorphic func- 
tion Z{\), 

(2.20) Z(X) = Y(\)e (e(x ^ 9 ^ a3 = ^=(a 3 + a 1 )\-i r7 ^{X)e- 9W ' 73 , 

i) Z(X) — ► I as A — > 00; 

(2.21) it) Z+(A) = Z_(A)G(A), G(A) = e<" T3 S{\)e-<" T3 , A £ j k , 

Z + (\) = o\Z-(X)<j\, X £ p. 

If 5(A) = I + sa± then G(A) = I + se ±2g a±. Our next goal is to transform the 
jump contour 7 to the contour of the steepest descent for the matrix G(A) — /. We 
denote by 7+ the level line Im(?(A) = const passing through the stationary phase 
point A = Ao = \J e~ iw x/6 and asymptotic to the rays arg A = ±ip This is the 
steepest descent line for e 2g . Let 7_ = Uj£j U a be the union of the level lines lj, 
j = 0, 1, 2, Im 5(A) = const, and a, Re g(X) = const, all emanating from the critical 
point A = — 2Ao. Among them, the level line l\ approaching the ray arg A = ^ 
(if argx = 7T, the level line l\ is the segment [— 2Ao,Ao]) is the steepest descent 
line for e 2g , while the level lines £0 and £2 approaching the rays arg A = — 4f and 
arg A = t?, respectively, are the steepest descent lines for e~ 2s . The level line a, 
Keg(X) = const, approaches the ray arg A = ir. 

Since the Stokes matrix 5i can be factorized, 

it is convenient to consider the following equivalent RH problems for ^(A): 

for arg a; £ [^,7r], the jump contour is the union of 7+ oriented from up to down 
and 7_ whose components are oriented toward infinity, see Figure |21 The jump 



Quasi-linear Stokes phenomenon for PI 



7 




Figure 2. A RH problem graph for s = 0. 



matrices are as follows: 

AG 7+ : 5(A) 

Aeli: 5(A) 

A€l Ul 2 : 5(A) 

(2.22) Aecr: 5(A) 

Remark 2.2. The jump contour for the RH problem l|2.22|) is decomposed into the 
disjoint union of the line 7+ and the graph 7 _, see Figure |5J For the boundary 
value argrr = it — 0, the level line l\ emanating from A = — 2Ao passes through the 
stationary phase point A = Ao and partially merges with the upper half of the level 
line 7+. To construct the RH problem however, it is not necessary to transform the 
jump contour precisely to the steepest descent graph. It is enough to ensure that 
the jump matrices approach the unit matrix uniformly with respect to A and fast 
enough as x — > 00 or A — * 00. 

As arga; £ [^,7r], introduce the reduced RH problem (s = s-i = 0) for the 
piece-wise holomorphic function 3>(A) discontinuous across 7_ only: 

i) lim A 1/2 (4=(<7 3 + ( 7 1 )A-3 ,T3 $(A)e- e,T3 - I) is diagonal, 

(2.23) ii) $+(A) = $_(A)5(A), A e 7- = U f =0 ,i, 2 <j U u. 
The jump matrix 5(A) here is dehned in l|2.22Jl . 



= 5_ 1 = 

= S+ = 

= 5_ = 
= —ia±. 







(0 
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Theorem 2.1. If argx G [^ ! "; 7r ] an d \ x \ * s lorge enough, then there exists a unique 
solution of the RH problem The Painleve function ya(x) corresponding to 

s = s_i = has the asymptotics Uq(x) — U 'e~ i7r x / '6 + 0{x~ 2 ) as x — > oo m i/ie 
above sector. 

Proof. Uniqueness. Since det5(A) = 1, we have det<J>+ = det$_, and hence 
dct $(A) is an entire function. Furthermore, because of normalization of <fr(A) at 
infinity, det $(A) = —1. Let $ and $ be two solutions of i|2.23[l . Taking into account 
the cyclic relation in l|2.8[l which implies the continuity of the RH problem for $(A) 
at A = — 2Ao, the "ratio" x(A) = $(A)4> _1 (A) is an entire function of A. Using the 
Liouville theorem and normalization of <E> and $ at infinity, we find x(A) = /, i.e. 

$(A) = $(A). 

Existence. Introduce an auxiliary function 

'v[(z) v' 2 (zT 
vi(z) v 2 (z) 



(2.24) = 



where the prime means differentiation w.r.t. z and 

(2.25) v x {z) = V2^e" /6 Ai(e l27r/3 z), v 2 (z) = -V2^M(z), 

with Ai(z) standing for the classical Airy function which can be defined using the 
Taylor expansion 0H1 EH , 

{ } [ ) ^ /3 nDto r (i)( 3fc ) ! 3l/3r (i)to r(|)(3fc + i)i • 

Asymptotics at infinity of this function and its derivative are as follows, 
(2.27) 

M(z) = ^ Z ^e-^ [^ lTi -^^ 

n—0 ^ 2 / ^ / 

Ai'(z) = ^^-^'{f; (-l)"3- 2 "(3n + |)g^,-^+C)(,-^)/ 2 )}, 

as z — > oo, arg z G (— 7r, 7r). 
It is worth to note that the function $o( z ) satisfies the linear differential equation 

(2.28) ^- = {za + + a-}i> . 

Using the properties of the Airy functions, we find that the products 

(2.29) = &o(z)S-, * 2 {z) = &i(z)S+, $ 3 (z) = $ 2 (z)S-, 
S± = I + io~± , have the asymptotic expansion 

(2.30) $ fe (z) = zi° 3 ±(a 3 + a^V^e^ 2 " 3 , 
as \z\ — > oo, argz e (— 7r + ^Lfc, -| + ^f-fc)), where 

v y w 4-; 2r(|)(2n)! V 6 " / 
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Let 7- = ctU.; = o,i,2 £j be the union of the rays £j — {z £ C: argz = =§-{j ' — 1)}, 
j = 0, 1,2, and a — {z G C: argz = 7r} all oriented toward infinity. This graph 
divides the complex z-plane into four regions: Q)q which is the sector between <x 
and £q, the sectors Uk, k = 1,2, between the rays £k-i and and the sector 0)3 
between the rays £2 and a. 




Figure 3. The model RH problem graph. 
Define a piece- wise holomorphic function $(z), 



(2.32) 
By cor 
(2.33) 



*k(«). 



By construction, this function solves the following RH problem, see Figure 03 



= I 



(2.34) 



*) ^(^^K^^K 
u) z<E7-: $+(z) = (0)5(2), 

zeii: §(z)=S+, ze£ U£ 2 : 

z £ a: S(z) = —itj\. 



0(z~ 
S(z) = S- 



Therefore the function <&(z) has precisely the jump properties of the function 
$(A). To find ^(A) with the correct asymptotic behavior at infinity, let us use the 
mapping 



(2.35) § z 3 / 2 = g(X) = f (A + 2A ) 5/2 - 4A (A + 2A ) 3 / 2 , 



z(A) = (-6A ) 2 / 3 (A + 2A )(1 - g^(A + 2A )) 2/3 , A = y/e~™x/6 



Within the disk |A + 2A | < R < 3|A | 



j 1 / 2 , the mapping l|2.35|l yields a holo- 



morphic change of the independent variable. Introduce a piece-wise holomorphic 
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(2.36) §(A) 



function $(A), 

'B(A)4(«(A)), |A + 2A |<i?, 
i (A + 2A )^ 3 ^ f (a 3 + a 1 ) e f( A ) CT3 , |A + 2A | > J?, 

S(A) = (-6Ao)-^(l-^)-* CT3 , 

where (A + 2Ao) 1 ^ 4 is denned on the plane cut along the level line a asymptotic to 
the ray argA = tt. Note that B(X) is holomorphic in the interior of the above disk 
|A + 2Ao| < R < 3 1 Ao | and thus does not affect the jump properties of <3E>(z(A)). We 
look for the solution of the RH problem l|2.23|l in the form of the product 

(2.37) 4>(A) = (7+ (4A 3 , - J£) C r+) X (A)$(A). 

Consider the RH problem for the correction function x(A). By construction, it is 
a piece- wise holomorphic function discontinuous across the clockwise oriented circle 
C of the radius R centered at — 2Ao and across the part of 7_ located outside the 
above circle (in fact, x(A) is continuous across a, see (|2.39|l below). The latter is 
divided by 7^ in four arcs: Cq between a and £$, C k , k — 1,2, between £ k -i and 
£k, and £3 between £2 and er, see Figure 0] To simplify our notations, let us put 

(2.38) A = A + 2A . 
Then the RH problem for x(A) is as follows: 

i) XW -> I, A ^ 00; 

U) X + (A)=x"(A)^(A), \€£, where 

A G £1, |A| > R: g(X) =1+ \<? 9 {o 3 - X 1 ' 2 o + + A _1 / 2 cr_), 

Ae4U4 \X\>R: G{X) =/+ie- 2 f(a 3 + A 1 / 2 ( 7 + -A- 1 / 2 ( 7_), 

A € u, |A| > R: G(X)=I, 

\X\=R, XeC k : G(X)^B(X)^ k (z(X))e^^(a 3 + a 1 )X-i^, 

(2.39) fc = Q,l,2,3. 

Taking into account the equations (|2.34|) , it is easy to check the continuity of the 
RH problem at the node points. Observing that, on the circle |A| = R = clxj 1 / 2 , 
< c < ^/3/2, we have z(X) — 0(|a;| 5 / 6 ) is large, we immediately see that 

(2.40) ||g(A) ^clAI 1 / 2 ^ 2 / 3 ) 172 !*! 172 ^ 372 , A e 4, k = 0,1,2, \X\ > R, 

where the precise value of the positive constant c is not important for us. Taking 
into account that, by the above reason, on the circle |A| = R, we may use for 
its asymptotics H2.30|l . the jump matrix G(X) has the asymptotic expansion 

(2.41) G{X) -1 = 



n=l 



^ 2T(i)(2n)! \l=kp_(i-Q n )\-l/z I±(_IT (1 _ 6n) 



= (-6A ) 2 / 3 (l- ^-) 2/3 A, A = A + 2A , \X\=R, 
5A ' 



Therefore, we have the estimate 

(2.42) \\g(X)-I\\<cR- 2 = c'\. 
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Figure 4. A RH problem graph for the correction function x(A). 

where the precise value of the positive constants c, d is not important for us. 

Now, the solubility of the RH problem l|2.39[l and therefore of l|2.23[l for large 
enough \x\ is straightforward. Indeed, consider the equivalent system of the non- 
homogeneous singular integral equations for the limiting value x + (A), i.e. 

(2.43, X+M -I- -L / *+«)(g"(C) -V 

or, in the symbolic form, \ + = I + K\ + . Here A + means the left limit of A on £ 
(recall, that the circle |A| = R is clock- wise oriented), and K is the composition of 
the operator of the right multiplication in Q~ 1 {\) — I and of the Cauchy operator 
C+. An equivalent singular integral equation for ip + ~ ^+ _ / differs from (|2.43|) 
in the inhomogeneous term only, 

(2.44) ifj + = KI + K^+. 

Consider the integral equation l|2.44|l in the space L2(t). Since t/" 1 (A) — / is small in 
L2(£) for large enough \x\, and C+ is bounded in L2(£), then ||JT||i 2 (^) < c|x| -1 / 2 
with some positive constant c, thus K is contracting and / — K is invertible in 
L2(£) for large enough |a;|. Because KI S L2(£), equation l|2.44|l for t/j + is solvable 
in L2(£), and the solution x(A) of the RH problem l|2.39|l is determined by ip + (X) 
using the equation \ + = I + KI + Kip + . 

Let us find the asymptotics of the Painleve function. Using 12.13fl and definition 
of $(A) (|2.36|l . the asymptotics of x(A) as A — > oo in terms of y and IK is as follows, 

(2.45) X (A) = /+ i(y - A - (4Ag - X) 2 )a 3 + j(4X 3 - M)cx_ + 

+ { 0{\- 2 ) 0(A- 3 / 2 ; 

On the other hand, in accord with the said above, the function x + is given by 
the converging iterative series, x + = S^Lo K n I. To compute the term K n I, we 
observe that the contribution of the infinite branches £k is exponentially small in x 
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due to estimate H2.40|) . Using the expansion (|2.41(l , we reduce the evaluation of the 
integral along the circle |A| = R to the residue theorem. Omitting this elementary 
computation, we present the final result: For large enough \x\, arga; € [%^,7r], the 
asymptotics of x(A) as A — * oo is given by 



(2.46) X(A)=/ + 



^(A^A- 1 ) 0(\q 4 \~ 1 ) / 

Comparing entries X2i(A) in l|2.45|) and l|2.46[) . we see that the Hamiltonian 
function !K = ^Cq(x) corresponding to the Stokes multipliers sq = s_i = is given 

by 

(2.47) K = K (x) = 4A 3 + C(Aq 2 ) = 4(-s/6) 3/2 + ©(x" 1 ). 

Next comparing entries Xn(A) in (|2.45ll and l|2.4(j[l and using 12. 4711 . we find the 
asymptotics of the Painleve function y = yo(x), 



(2.48) y = y (x) = A + (D(\- 4 ) = ^x~/6 + G(^ 2 ). 

Recall that Ao = (er n 'x/6) 1 ^ 2 where the main branch of the root is taken. □ 

Let us go to the case of the nontrivial s_i described by the RH problem 12.22JI . 
We look for the solution ^>(X) in the form of the product 

(2.49) *(A) = (I-(K- M )a + )X(A)$(A), 

where $(A) is the solution of the reduced RH problem (|2.23|1 and !Ho (j2.47|) is 
the Hamiltonian function l|2.14(l corresponding to the Painleve transcendent yo(x) 
(|2.48ll . Using 1)2. we find the asymptotics of AT (A) as A — * oo, 

(2.50) X(X) = (7+(5C- JfoW)** -1 = 

~ 2A ^ ~~ ' ~ ~~ X ~ Ji o) (J - + I o(X~ 2 ) 0(A~ 3/2 ) 

Thus we arrive at the RH problem for the correction function X(X) on the 
steepest descent line 7+, 

i) X(A)->J, X —> 00, 

(2.51) it) X+(A)=X_(A)S(A), Ae 7+ , 

9(A) = $(A)5_i$- 1 (A). 

Note, $(A) is continuous across 7+ and therefore holomorphic in some neighborhood 
of 7+ as argx £ [^-,n], \x\ is large enough. 

The jump matrix on 7+ can be estimated as follows, 

(2.52) ||g(A) - 7|| < c | a _ 1 | e -i2»/« 3 V«| !B |V« co(l(4(arg!e _, ) ) e _ 2 s/V/*k|V«|A-Ao| a < 

Here c is some positive constant whose precise value is not important for us, and 
Ao = (e'^a;) 1 / 2 is the stationary phase point for exp(<?(A)), see Q2.18jl . Estimate 
(|2.52[1 yields the estimate for the norm of the singular integral operator JC in the 
equivalent system of singular integral equations, Af_ = I + JCX-, 

(2.53) ||/C|U 2(7+) < c/]a_ x |e-* a " / * 3l/ *l»l ,,/ * ct -(t<~« SB — d > 0. 

If |x| is large enough and argx € [nf + £ ; 7r], e > 0, then the operator JC is contracting 
and the system A"_ = I + JCX- is solvable by iterations in L 2 (7+), i-e. X- = 
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Y^Lo fc n X-. However, to incorporate the oscillating direction argx = 4p in the 
general scheme, we use some more refined procedure. 

Theorem 2.2. If sq = 0, argx £ [^-,tt] and \x\ is large enough, then there exists 
a unique solution of the RH problem \2. 1(^1 - 12.12^ . The corresponding Painleve 
function has the asymptotics 
(2.54) 

V(x) = yo(x) + ^2- 11 /8 3 - 1 /8( e -i 7 r a; )-l/8 e -|2"/43V4 (e -- :c )V4 + 

where yo(x) ~ \J e _i7r a:/6 is the solution of the Painleve equation for sq = s_i = 0, 
si = s 2 = S-2 = i- 

Proof. It is enough to prove the solubility of the RH problem l|2.51|l . 

Using for $(A) the expressions (12.37(1 with (|2.36() and the estimate (|2.46|) to- 
gether, we find the asymptotics of the jump matrix S(A), 

(2 551 qfAW+ia ie 2 *^ 1 + (9(A ° 2)r ~ /2) -^ 1/2 + ^(V 2 ) 

(2.55) b{\)-l + 2 s- ie [ x -i/2 + 0{x -2 x - 1) _ 1 + 0(A -2 A -i/ 2) 

AG7+, A = A + 2A . 

Consider the following model RH problem, 
i) P(A)->7, A^oo, 

(2.56) ii) P+(X) = P_(A)§(A), Ae 7 +, 

1 -(SAo) 1 ^ 

This problem is solvable by the following quadrature, 

(«-, p W =, + k iA £^c( (3Ao ;_ 1/2 - (3 Y' 2 

We look for the solution X(X) of the RH problem ((2.51(1 in the form of the 
product, 

(2.58) X(\) = Q(\)P{\). 
The correction function Q(X) satisfies the RH problem 

i) Q(A) — ► J, A — > oo, 

ii) Q + (X)=Q.(X)W(X), Ae 7+ , 

(2.59) W{\) = P_(A)S(A)S(A)- 1 Pr 1 (A). 

Using 1(2.55(1 - 1(2.57(1 . we find the estimate for the jump matrix W(A) on 7+, 

(2.60) VF(A) =/ + 0(s_ ie 29 (A- A )Ao 1/2 ), A e 7+, 

Our next steps are similar to presented in the proof of Theorem l2.ll Consider the 
system of the singular integral equations for Q + (\) equivalent to the RH problem 
(12.59(1 . Q + = I + 3CQ+. Here the singular integral operator % is the superposition 
of the multiplication operator in W — I and of the Cauchy operator C+ . Because 
the Cauchy operator is bounded in £2(7+), the singular integral operator % for 
large enough argee £ [^,7r], satisfies the estimate 

(2.61) ||3C|U 2(7+) < c | J! _ 1 || a: |-V2 e -i2 I ^3V«|«|V« COB( i (arga! _ T))j 



S(A) = / + K ie 29 V(3Aorl/2 _ x 
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with some positive constant c whose precise value is not important for us. Thus 
equation = XI + X(+ for the difference £+ := Q+ — 7 is solvable by iterations in 
the space £2(7+) for large enough \x\. Solution of the RH problem (|2.59|) is given 
by the integral Q = I + XI + X( + . This implies the asymptotics of Q(A) as A — > 00, 



(2.62) Q(\)=I+J-J (i + OiXIiOfiil-W- 1 ^ 



= I + O^s-rx- 1 ' 2 exp(-|2 11 / 4 3 1 / 4 Re (e^) 5 / 4 )) . 

Now let us find the asymptotics of the Painleve function y(x). Using Q2.58p . 
(|2.57|l and the estimate (|2.62|) , we find 

(2.63) 

X(\) = I + i_L 2 - 19 / 8 3- 1 / 8 (e— zyl/te-i* 1 '**'*^*'*)*'* (/ + 0(x^)) x 

1 -2- 1 /4 3 l/4 (e -™ a;) l/4N 
2 l/4 3 -l/4( e -i Wa .)-l/4 _ x 

Comparing (|2.63|) and (|2.5()|l , we conclude that the Hamiltonian function for so = 

is as follows, 

(2.64) 

X(x) = X (x) - f-i2- 17 / 8 3" 3 / 8 (e'"x)- 3 / 8 e-^ 2ll/43l/4 ( e '" a; ) 5/4 (l + 0(x- 1/8 )), 

while the Painleve function is given by 
(2.65) 

y 



{x) = yo{x) + £^ 2 -ll/8 3 -l/8 (e -^ :c) -l/8 e -l2"/43 1 /4 (e - a:) 5/4 (i + 0(a .-3/ 8))) 



where 5Co(^) and yo^) & re the Hamiltonian and the Painleve functions, respectively, 
corresponding to sq = ,s_i =0. □ 

2.2. Other degenerate Painleve functions. Applying the symmetry (|2.16a|l to 

the solution l|2.54|l and changing the argument of x in 2n, we obtain 

Theorem 2.3. if so = and \x\ — > 00, argx £ [n, t?-], then the asymptotics of the 

Painleve first transcendent is given by 

(2.66) 

y(x) = y x {x) - ^-^Z-V^e-^xy^e-^'^^^il + O^)), 

where yi(x) ~ yj e~ i7T x/6 is the solution of the Painleve equation for sq = Si = 0, 
s-i = s 2 = S-2 = i. 



The solutions yo(x) and yi(x) = yo(e 2m x) are meromorphic functions of x £ C 
and thus can be continued beyond the sectors indicated in Theorems 12 . 21 and 12 . 31 
To find the asymptotics of yi(x) in the interior of the sector arga; £ [^-,tt], we 
apply (|2.54|) . Similarly, we find the asymptotics of the solution yo(x) in the interior 
of the sector argee £ [it, ^] using l|2.66|l . Either expression implies 

Corollary 2.4. If \x\ -> 00, argx £ [^, 7 -f], then 
(2.67) 

Vl{x) - y (x) = -4 = 2- 11 /8 3 - 1 / 8 (e-^ 2 ;)- 1 /8 e -|2 11 / 4 3 1 / 4 (e — + 0(x -3/ 8)) 

\/7T 
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Applying symmetries (|2.16(l to yk{x), k = 0,1, we find the solutions yk{x) cor- 
responding to the Stokes multipliers Sk — Sk-i = 0, 

V2n{x) = e l ^- n y (e l2 S- n x) for s 2 „ = s 2n -i = 0, 

(2.68) 2/2n+i(z) =e l ^ n yi (e l ^ n x) for s 2 „ +1 = s 2n = 0. 

Since there is one-to-one correspondence between the points of the monodromy 
surface and the Painleve functions, the first of the identities (|2"lH]l . Sk+5 = Sk, 
implies that y n+5 (x) = y n (x). 
Using Theorems O and E31 

we find that 

(2.69) y 4n (x) = y An+5 {x) = je-™% + 0(x~% 

\x\ > oo, arga; € [§ - ^n,7r - ^n], 

(2.70) y in -2(x) = 2/4„+ 3 ^) - + 0(^ 2 ), 

\x\->oo, &rgx£[f-fn,^-fn}. 
The symmetry ()2.16jl with the definition Ij2.68j) applied to Ij2.67|l yields 
Corollary 2.5. // \x\ oo and argz € - ^n, ^ - ^n], i/ien 

(2.71) ^n+iO) -V2n{x) = 

= e i i+^" 2 _ 11/83 _ 1/8 i7r+i% 1/8c _i 2 n/4 3 i/4 (e -i. + ^n x) 5/4 ;; 

x (l + 0(x~V 8 )), 

On the one hand, equations l|2.67|l . (|2.71|l constitute the quasi-linear Stokes 
phenomenon for the Painleve first equation. On the other hand, these equations 
give the asymptotic description of the degenerate Painleve functions beyond the 
sectors in l|2.69l) and H2.70(l . Observing that the difference l|2.71|l is exponentially 
small in the interior of the indicated sector, we conclude that the asymptotics (|2.69(1 
and H2.70JI as \x\ — > oo continue to wider open sectors, 

(2.72) y in (x) = yV^f + 0{x~ 2 ), argz e (e - § - f n, If - % n - e), 

(2.73) y in - 2 (x) = -y/e-™ f + 0(x- 2 ), aigx e (e + f - f n, ^ - f n - e), 
where e > is an arbitrary small constant. 

Remark 2.3. The solutions y n {x) (|2.68|) corresponding to the trivial values of two 
Stokes multipliers s n = s n _i = are the most degenerate among the Painleve 
transcendent since they behave algebraically in four of five sectors arga; e (-§ + 
f^f + ?fk), k = 0,±1,±2, see $FT3 . l|2~75| l. Nevertheless these solutions are 
transcendent, since their asymptotics as \x\ — > oo within the remaining fifth sector 
involves the elliptic function of Weierstrafi, look for more details in [30]. Moreover, 
the fact that the asymptotics of y n (x) is not elliptic in four sectors uniquely de- 
termines the values of all the Stokes multipliers Sk- Thus the asymptotics Ij2.72|l . 
(|2.73Jl uniquely determine the degenerate solutions y n (x). 

Remark 2.4. The asymptotics of less degenerate solutions corresponding to s n = 
and ,s„ + i + s n -i = i can be found applying the symmetries Ij2.16b() to equations 
(l2~5l> and iglSJ). 
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3. Coefficient asymptotics 

Using the steepest descent approach, cf. 02], we can show the existence of the 
asymptotic expansion of y n (x), n G Z, in the negative degrees of x 1 / 2 . Further 
elementary investigation of the recursion relation for the coefficients of the series 
allows us to claim that the asymptotic expansion for y n (x) in (|2.72|) . 1|2.7L)[I has the 
following form: 



(3.1) y f (x) = a(-f ) 1/2 £ a k a k (-xr 5 ^ + 0(x 

k=0 



= a(-l) 1/2 ^2^ k (-x)- 5k + ^(-x)- 2 Y / a 2 k+ i(~x)- 5k + 0(x-n, o* = l, 

fc=0 k=0 

where coefficients a k are determined uniquely by the recurrence relation 

fe 

(3.2) a Q = 1, afe + i = 25k ^ 1 a,k — \ ^ a m ak+i- m - 

m—l 

Several initial terms of the expansion are given by 

(o Q\ — rrJ Wfi/l I 49 4412401 , 245229441961 , n(„-20\\ 

(6.6) y } (x) - o^-x/'oy + ^ggpr - 1179648:c iu + 100663296^ + <~>(x jj>- 

1_ /i _ 1225 r 73560025 _ 7759635184525 , -20\ \ 

48a: 2 \ 192a; 5 "l" 49152a; 10 3538944a; 15 If- 

Our next goal is to determine the asymptotics of the coefficients a& in l|3.1|l as 
k — ► oo. With this purpose, let us construct a sectorial analytic function y(t), 

(3.4) argi e [-f (n + l),-f n ] . y (t) = y 4n (e l7r t 2 ), n = -2, -1,0,1,2. 

The function t/(t) has a finite number of poles all contained in a circle |i| < p and 
is characterized by the uniform asymptotic expansion near infinity, 

, oo 

(3.5) y(t) = _^ afct -^ + 0(r~). 
Let y^ N \t) be a partial sum 

(3.6) yW(t) = ±Y,a k tS\ 

V6 ^o 

and y( N \t) be a product 

oo 

(3.7) = i 5Ar - 2 V6(y(t) - yW(i)) = f" 1 £ a fe+A rt- 5fe + £)(r»). 

fc=0 

Because t 5N ~ 2 y( N \t) is polynomial, the integral of v^ N '(t) along the counter-clock- 
wise oriented circle of the radius \t\ = p satisfies the estimate 
(3.8) 

/ v {N \t)dt <p 5N - 2 V&<f> \y{t)\dl < V62irp 5N - 1 m&x\y(t)\ = Cp 5N 
l\t\= P J\A=p l*l =p 

with some positive constant C whose precise value is not important for us. 
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On the other hand, inflating the sectorial arcs of the circle |t| — p, we find that 

2 

(3.9) / v^ N \t)dt=l v^(t)dt+f] [ (vf\t)-v [N \t))dt. 

J\t\=p J\t\=R n =-2 Jel2i ~ n (P' R ) 

Because v^ N \t) = t _1 ajv + 0(t~ 6 ), the first of the integrals in the r.h.s. of l|3.9ll is 
computed as follows: 

(3.10) I v (N) {t)dt = 2Tria N + <D(R- 5 ). 

J\t\=R 

Remaining integrals in l|3.9|) are computed using definitions <|3.4[l 13. 7[) and (|2.68() 
with the identification y_4(x) = yi(x) and the formula l|2.67[l together, 

(3.11) 

2 

£ J (vi N \t)~v W (t))dt = 5V6 J t 5N -Z(y^(e™t 2 )-y (eS))dt = 

= z ^ 2 - 11 /8 3 -l/8 f ^-f e -*2" / *3VV/' (1 + 0(t -3/4 ))dt = 

v n J 

(p,R) 

= 2i^S-(\2 11 ' i ?, 1 / A y 2N T(2N - i)(l + 0(N- 3 / 10 )) + 0{p m ~i)+ 

+ <D(e-^ 11/i3l/iR5/2 R 5N ^). 
Thus, letting R — oo, we find the asymptotics of the coefficient in l|3.1|l as 

N -> GO, 

(3.12) 

a N = -^^(p^/y 2N T(2N - |)(1 + 0(7V- 3 A°)) + 0( p ^), N ^ 

Remark 3.1. The presented asymptotic formula shows a remarkable accuracy: ne- 
glecting in H3.12f> error terms, we find an approximation to ajv with the relative 
error not exceeding 2% for N = 4 and 1% for N = 7. Furthermore, for the initial 
set of N = 1, 2, . . . , 7, the relative error decreases approximately as iV _1 which is 
significantly better than estimated. 
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